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Abstract 

The interaction between singular and regular fields is considered for Lorentz- 
invariant scalar and vector wave equations. The singular field is generated by a Dirac 
source term. Its dynamics are deduced from the total field Lagrangian. At non- 
relativistic speeds, the resulting equations of motion are those of a mass in a scalar 
potential. Using this method we deduce the relationship between source amplitude 
(scalar gravitational mass) and dynamic mass (inertial mass). Generalising this 
method implies Lorentz forces for charge singularities in the electromagnetic field 
and describes the dynamics and interaction of hypothetical magnetic monopoles. 

1 Introduction 

The method presented builds on the action principle which offers a general tool to 
describe both particle and field dynamics. This principle represents one of the great 
unification in physics and applies equally well to Newtonian mechanics pQ, general 
relativity and electromagnetic fields Using this single principle it is possible 
to deduce the Euler-Lagrange equations describing the evolution of one system or 
of many interacting systems. However describing these interactions requires the 
assumption of an additional action term. In this paper, we show that it is possible 
to eliminate this additional assumption and deduce the reaction dynamics of the 
singularity. In the case of scalar field singularities we derive Newton's equation of 
motion together with the equivalence between Newtonian gravitational mass and 
inertial mass. For charges and magnetic monopoles we deduce their associated 
Lorentz force. 

Singularities, vortices and topological defects have received considerable atten- 
tion over the years; they determine the properties of many interesting materials such 
as vortices in superconductors superfluids [I] and two-dimensional condensate 

dislocations and defects in solid [Hj and liquid crystals |7j, optical phase singu- 
larities and optical vortices [B] and even topological defect in cosmology [§]. In this 
paper, we use singularities together with the field action principle to deduce their 
dynamical behaviour. The method presented here is general and can be applied to 
all of the fields mentioned above. 
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For clarity reasons we treat first the simplest case possible that being a Dirac 
source singularity in the Lorentz-invariant scalar wave equation. This leads to equa- 
tions of motion similar to Newtonian mechanics (section 2). The method is gener- 
alised for other singularities and other fields such as a charges (section 3) where the 
electromagnetic field implies the Lorentz force. In the last section, we apply this 
principle to deduce the equations of motion of hypothetical magnetic monopoles 
where the dynamics concurs with the one obtained through charge-monopole sym- 
metry considerations |10|. 



2 Scalar field singularity: Inertia 

Before considering the anatomy of the field singularity, let us start by deducing the 
Euler-Lagrange equation for a scalar mass-less field U described by the following 
action 

Su = J C{U) d 3 x dt =\J ^( d tU) 2 - (VC/) 2 d 3 x dt (1) 

where C(U) corresponds to the Lagrangian density, V = (d x ,d y ,d z ) to the nabla 
operator and c to the speed of light. The action principle implies that the scalar 
field U is solution of the standard homogenous wave equation 

V 2 U - \d 2 t U = (2) 



c 



where V 2 stands for the Laplacian operator. The next step in the singularity dynam- 
ics theory is to define a singular source term that maintains the Lorentz-invariancc 
of the supporting field equation. Starting from the wave equation (0, we can define 
the stationary singular field U s as the time independent solution in the presence of 
a Dirac source at the origin of the coordinate system 



V 2 U S - \d 2 t U s = m s 5(r) (3) 



where m s corresponds for the source amplitude. To define the Lorentz-invariant 
singularity, we consider equation in a reference frame having the constant velocity 
—v with respect to the stationary singularity. In this reference frame the singularity 
has the velocity v. The transformation between the two reference frames is 1 1 1 j : 

r'|| = -y (r || + vt) 

r'± = r x (4) 
t' = l{t+—. 



where 7 = (1 — v 2 /c 2 ) -1 / 2 and the prime stands for the coordinates in the moving 
reference frame. The subscripts || and _L represent the vectorial components parallel 
and perpendicular to the velocity v. Applying these transformations to equation Jjyl 
leads us to a constant velocity moving singularity. Beforehand, we remark that the 
singularity can be decomposed into two components, 8{r) = <5(ri|)5(r_|_). Further, the 
wave equation and the scalar field U s is invariant with respect to the transformation 
that is it remains unchanged in the transition from one reference frame to another. 
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It is only the singularity which changes in the moving reference system 



V' 2 U s -^d 2 U s = m s 5( 7 (r'|| -vt'))S(r' ± ) 



l-^r'n-vt'))^!) 

cr 



m s yi-^(r'-r' s ) 

where r' s = vt' is the position of the singularity in the moving reference system. We 
have used the scaling property of the single variable Dirac distribution |a|<5(aa;) = 
S(x). 

We can now write the equation of a moving singularity on an arbitrary path r s (£) 
with the velocity v = r s . The singular field U s is solution of 
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V 2 U S - ^d 2 U s = m s Vl-v 2 /c 2 *(r - r s ). (5) 

and corresponds to a retarded potential. 

In the following, we are considering the action principle in the case of the linear 
superposition of the singular field as defined in equation J5J and a given external field 
Uq solution of @. We remark here that the only free parameter of the system is the 
position of the Dirac singularity and its velocity. The Euler-Lagrange equation of 
motion is deduced when the action is stationary with respect path of the singularity. 
The total action integral may be expressed as 



S Uo+Us {r s (t)) = J d 3 xdtC(U Q + U s ) 

= \J d?x dt ^(0 t U ) 2 - (Vt/o) 2 + ^{dtU s ) 2 - (VU S ) 2 

+ J d 3 x dt ±(d t U )(d t U s ) - (VU )(VU S ). 

This action integral can be decomposed into two parts. The first integral has no 
cross-terms between the external field and the singular field. This integral is inde- 
pendent with respect to the position of the singularity and its variation is zero. This 
leaves only the second integral AS S for the action principle 



AS S = J d 3 x dt U V 2 U S - ^U d 2 U s 



where we have integrated the action by parts and assumed the singular field to be 
zero on the boundary of the integration region. Next, we use equation Q to replace 
the singular field by the Dirac singularity and finally we integrate in Lagrangian 
density in the spacial domain. 

AS S = J d 3 x dt U m s ^l - v 2 / 2 S(r - r s ) 

dt U m sy Jl-v 2 /c 2 
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The equations of movement can be obtained by using the conventional particle Euler- 
Lagrange differential equation where the position and velocity correspond to those 
of the singularity. Indeed we have 



d_ 

dt 



rriiV = -m g VU (6) 



where we have defined an effective inertial rrii and gravitational m g mass 

U m s 

rrii 



yi-v'/c 2 (7) 



m g = m s \/X — v 2 /c 2 . 

One can also define a differential relationship between the two masses by taking 
the scalar product between the velocity and equation ©. This leads to a energy 
balance equation that generalises the famous equations E = mc 2 to the dynamic 
domain. 

d 2 d TT 

-m lC = m g -U (8) 

Indeed, the energy balance equation shows that the variation of the inertial mass is 
proportional to its potential energy variation. 

Finally, we can deduce Newtons equation of motion for non-relativistic velocities 
and far from any other singularities 

±m s v = (9) 

where the scalar field gauge is chosen to be f7o(oo) = c 2 . 



3 Charge singularity: Lorentz force 

In order to deduce the laws of motion of a charge we have to consider a singularity of 
the electromagnetic field described by Maxwell's equations. We proceed in a similar 
way to the scalar singularity case. In a first instance, we deduce the homogenous 
Maxwell's equations using the action principle. In a second step, we introduce a 
singular source to these equations making sure that their symmetries are maintained 
(i.e. Lorentz-invariance) . Then, we write down the Lagrangian integral of the total 
field composed of an external field and the singular field. Finally, the action principle 
is applied to the total Lagrangian integral and the dynamics of the singularity is 
deduced. 

The electromagnetic field action is expressed as: 

Sera = J ( yE • E - J-B • d 3 X dt (10) 

where the electric and magnetic fields are defined by 

E = -VV-d t A 

B = VxA (11) 
= V-A + e of i d t V 
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and where V and A are the electric scalar and magnetic vector potential respectively. 
The vacuum permittivity and permeability is given by eo and (Iq. The third equation 
corresponds to Lorentz gauge. Applying the action principle with respect to the 
potentials lets us deduce Maxwell's equations in absence of charges: 



VxE-ftB = (12a) 

Mo 



-VxB-foftE - (12b) 



V ■ e E = (12c) 

V B = 0. (12d) 

Like in the scalar case, we consider a single stationary Dirac singularity qS(r) in 
equation of the divergence of the electrical field l|12c|l . The effect on the electromag- 
netic field of this source is equivalent to a stationary point charge q at the origin. 
The singular charge fields E q and B q are solutions of the in-homogenous Maxwell's 
equations: 

-V x E q - d t B q = 
— V x B q - e a t E q = 

Mo (13) 

V ■ e E q = qS(r) 

V • B q = 

The main difference with the scalar case is in that the potential fields V q and 
A q change in the transition from one reference frame to another. Indeed, using the 
moving frame defined by (J2J we have for the potentials: 

c 2 J 

(14) 



A q|| = 7 (^A q 
A' — A 

K=^+V-A q ) 

where the speed of light is given by 1 = c 2 eo/xo- Combining the transformations J5J 
and 114(1 together with the definition of the electric and magnetic fields ()ll(l we can 
deduce the transformation relationship between these fields 



E q|| 


= E q|| 


E qX 


= 7(E q± - 


B «u 


= B q|| 


B qX 


= 7 ( B q± 



XB q ) 



V X E q 



(15) 



This implies the following transformation for equations 113|) : 

-v'xe;-p; = o 

Mo 

V'-e E q = g 7 5( 7 (r'||-vt / ))5(r / ± ) 

V • B q = 
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Applying the scaling property of the single variable Dirac distribution |a|<5(ax) = 
d(x) and introducing the definition of the position of the singularity r' q = vt' in the 
moving reference system we get: 

v'xe; + 9,b; = o 

-V'xB;-e„p; = qvS(r'-r' q ) 
Mo 

V'-eoE^ = q6(r'-r' g ) 

V • B q = 

where the relativistic coefficient 7 in the source terms cancelled out. The reason for 
this is the covariant behaviour of the electromagnetic potentials as opposed to the 
invariant scalar potential in equation 

This procedure gives us the form of the singularity following an arbitrary path 
r q (t) with the velocity v = r q . The resulting vector singularity is equivalent to a 
charge and its associated current. The singular charge fields E g and B,j are solutions 
of the in-homogenous Maxwell's equations: 

- V x E q - <9 t B q 

— V x B q - 6 d t E q 
Mo 

V • e E q 
V B q 

To deduce the equations of motion of this moving charge singularity we define 
the total field as the linear superposition of the singular fields E q and B q and the 
external fields Eo and Bo: 

E = E q + Eo 
B = B q + Bo 

where the external fields are decomposed onto the external scalar and vector poten- 
tials 

E = -W - d t A 
B = V x A 
= V ■ A + 6oMo<9 t V . 

Introducing the total fields E and B into the electromagnetic Lagrangian 
infers the following singularity action: 

AS sq = J ^e E ■ E q - -i-B ■ B q ^ d 3 x dt (17) 

where we have kept only the terms that vary with respect to the trajectory of the 
singularity. Indeed, the terms in Ej • Ei and Bi • B{, where (i = q, 0), are not 
dependent on the position of the charge singularity. Rewriting the external field 
with the help of the scalar Vq and vector potential A allows us to integrate by 







= qvS(r - r q ) 



qS(r - r q ) 
0. 



(16) 
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parts: 



AS sq = j (e (-Wo - d t A ) ■ E q - i- (V X A ) • Bq) cPx eft 

= y fe (F V • E q + A • 9*E q ) - ±- (V x B q ) • A j rf 3 x A 

where we have assumed the fields E q and B q to be zero on the boundary of the 
integration domain. Using equations (|16f) we can further simplify the action of the 
singularity 

AS sq = J (qV S(r - r q ) - qA ■ vS(r - r q )) d 3 x dt 

= J (qV - qA ■ v) dt 

which leads directly to Lorentz forces for a charge in the external field defined by Eo 
and Bo- To convince ourselves of this we can rigidly link a scalar mass singularity 
introduced in previous section to a charge singularity. The total action of this double 
singularity is: 



AS S + AS sq = I ( U m s \/ 1 - J + qV - qA -vjdt (18) 

which lead to the following equation of movement 

^m,v = -m g VU + qE + qv x B (19) 
where n%i and m g are the defined by JJJ). The energy balance equation is 

^m lC 2 = m g -^[/ + gv-E . (20) 

4 Magnetic Monopole singularity 

Magnetic monopoles are hypothetical particles that behave in a similarly to charged 
particles with magnetic and electric fields inverted (|1()|). The monopoles correspond 
to sources of the magnetic field in equation (|12dl) . In the following we are deducing 
the equations of motion of magnetic source singularities by considering a Dirac 
source q m S(r) in equations (|12d|) where q m corresponds to the monopole strength. 

Before proceeding further, we have to generalise the potential decomposition of 
the electromagnetic field. Indeed, the potentials that we used for the charged singu- 
larities does not allow the presence of sources of the magnetic field. Conventionally 
the magnetic field is by definition a rotational field. Therefore, we define a symmet- 
ric electromagnetic potential decomposition that includes an additional magnetic 
scalar and electric vector potential. 

E = - W 9 - d t A q + —V x A m 

B = V x A q + fi Q VV m + /! <9 t A m (21) 
= V • A q + eofiodtVq 
= V ■ A m + e Hod t V m 
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where the subscripts (q, m) stand respectively for the charge and monopole poten- 
tials or fields. Using all the newly defined potential definitions we can apply again 
the action principle the electromagnetic action defined by equation l|10|) . The re- 
sulting evolution equations are identical to equations i|12|) . The only difference is 
the possibility to define a monopole source density in equation l|22d|) . 



- V x E - d t B = (22a) 

-VxB-e 9 ( E = (22b) 
Mo 

V-e E = (22c) 

V B = 0. (22d) 

For the monopole singularity, we consider a single stationary Dirac singularity 
at the origin. The singularity constitutes the free term of equation i|22d(l and acts 
as a source of magnetic field. The singular fields E m and B m are solution of: 

-V x E m - 8 t B m = 

— V x B m - e d t E m = 

Mo (23) 
V • e E m = 
V B m = q m 6(r) 

In this case there are four potential fields V sq , V sm , A sq and A sm and the 
electric and magnetic field is defined by equations i|21|l . Taking the same frame 
transformation defined by J2J we have for the potentials: 

A ' - y(A + ^) 
A sq|| — 7l A sq|| -1- ^ ) 

A' — A 

^=7(K 9 +v-A sq ) ^ 
A' - ^CA „ + vVsm ) 

A' — A 

C = 7^m+v - A sm ). 

The implied relationship between the electromagnetic fields in the two frames is 
identical to the transformation for the charges only field. 

E m|| = E m|| 



E' m± = 7 (E m± - v x B m ) 
v x E 



B' m ||=B m || (25) 



Combining all the above transformations implies the singularity in equations (|23J) 
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to transform as: 

-V'xE' m -9 ( < = gro v7<K7(r'||-vi'))<5(r'±) 

-V'xB^- £o a t < = 

Mo 

V'-eoE^ = 

V'B' m = g rn7 J(7(r'|| - vOWrM 

Applying the scaling property of the single variable Dirac distribution |a|(5(aa;) = 
and introducing the definition of the position of the singularity r' m = vt' in 
the moving reference system we get: 

-V'xE^-3A = Qmv6(v' - r' m ) 

iv'xB' m -e ^E' m = 

Mo 

V'-eoE^ = 

V'B' m = g m( 5(r' - r' m ) 

where the relativistic coefficient 7 in the source terms cancelled out just like in the 
charge singularity case. 

Consequently, the electromagnetic singular fields E m and B m are solution of 

-e V x E m - — <9 t B m = q m v5(r - r m ) 
Mo 

—V x B m e 9 t E m = 

Mo (26) 
V • e E m = 

V B m = q m 5(r - r m ). 

Mo 

for a magnetic monopolc singularity on an arbitrary path r m (t) with a velocity 
v = f m . 

The equation of motion of the monopole singularity is deduced by defining the 
total field as the linear superposition of the singular fields E m and B m and the 
external fields E and B : 

E = E m + Eo 
B = B m + B 

where the external fields can be decomposed onto four potential fields V~o q , Vo m , A 0q 
and A 0m 



E = -W , - d t A 0q + —V x A 0m (27a) 

£o 

B = V x A 0q + Mo VVb,„ + Aio<9 t A 0m (27b) 

= V • A 0q + eoHod t Va q (27c) 

= V-A 0m + e Mo<9 t Vb m . (27d) 
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Introducing the total field in the electromagnetic Lagrangian (eq. I1(J|) infers the 
following singularity action: 

AS sm = J f e Eo • E m - -^B • B m ^j d 3 x dt (28) 

where we have kept only the terms that vary with respect to the trajectory of the 
singularity. Indeed, the terms in Ei ■ Ei and Bi ■ Bi, where (i = m, 0), are not 
dependent on the position of the charge singularity. Rewriting the external field 
with the help of the potentials defined in equations H27a|) and (|27b|) allows us to 
perform an integration by part 



AS sm = J (eo (- W g - <9 t A 0q + V x A 0m ) • E m 

— (V x A 0q + VV Qm + d t A 0m ) • B m ) d 3 x dt 
Mo / 

= f (eo (V 0q VE m + A 0q <9 t E m + A 0m V x E m ) 

-— (A 0q V x B m - V 0m VB m - A 0m a t B m ) J d 3 x dt 

Mo / 

where we have assumed the fields E m and B m to be zero at infinity i.e. the bound- 
ary of integration. Using equations l|26|) we can further simplify the action of the 
singularity 

AS sm = J (q m V 0m 6(r - r q ) - q m A 0m ■ vS(r - r q )) d 3 x dt 

= J (q m V 0m - q m A 0m ■ v) dt 

which leads directly to a Lorentz type forces. Like in the case of the charge singu- 
larity, we can rigidly link a scalar mass singularity introduced in previous section to 
a monopole singularity. The total action of the double singularity is: 



A5 S i- A«S„„ / ( U Q m s \/l - ^ + q m V 0m - q m A 0m ■ v j dt (29) 



— miV = -m g VUo + q m B 0m + q m v x E 0m (30) 



which lead to the following equation of movement 
d 
It 

where rrii and m g are the defined by J7J) and 

Eom = — V X Aom 

eo 

Bom = MoVVo, n + (lodt Aom 

correspond to the magnetic monopole part of the electromagnetic field. We remark 
here that the force acting on the monopole is only due to the monopole part of the 
electromagnetic field. Consequently, there is no direct interaction between charges 
and magnetic monopoles. 

The energy balance equation is 

d d 

-j^rriiC 2 = m g — U + q m v-B Om . (32) 
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5 Conclusion 



The field generated by a singularity gives rise to the laws of motion of the singularity. 
In other words, the action of the singularity on the field implies its dynamic reaction 
to the field. We have shown this to be the case for neutral bodies, charges and 
hypothetical magnetic monopoles. Using our method we have deduced not only the 
dynamics of masses, but also the relationship between the "gravitational" source 
strength and the associated inertial mass. The charge singularity was shown to be 
the origin of the Lorentz force. The magnetic monopoles have an equivalent Lorentz 
force. 

The conventional Lorentz force does not involve the magnetic monopole field. 
Reciprocally, the monopole Lorentz force does not involve fields originating from 
charged particles. Consequently, charges and magnetic monopoles do not interact, 
making the discovery of monopoles using charges difficult. A major question to be 
addressed concerns the role played by electromagnetic radiation. Can it provide the 
link between charges and magnetic monopoles? 
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